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Abstract 

The HyperCP collaboration has observed three events for the decay S + — > pfi + fi~ which may be in- 

• terpreted as a new particle of mass 214.3 MeV. However, existing data from kaon and S-meson decays 

severely constrain this interpretation, and it is nontrivial to construct a model consistent with all the data. 

, In this letter we show that the "HyperCP particle" can be identified with the light pseudoscalar Higgs 

boson in the next-to-minimal supersymmetric standard model, the A\. In this model there are regions of 

Q_i| parameter space where the A® can satisfy all the existing constraints from kaon and S-meson decays and 

O • mediate S + — > pfi + u~ at a level consistent with the HyperCP observation. 
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Three events for the decay mode E + — > pfi + fi~ with a dimuon invariant mass of 214.3 MeV have 
been recently observed by the HyperCP Collaboration [1] . It is possible to account for these events 
within the Standard Model (SM) when long-distance contributions are properly included • How- 
ever, the probability of having all three events at the same dimuon mass, given the SM predictions, 
is less than one percent. This suggests a new-particle interpretation for these events, for which the 
branching ratio is (3.llfJ ± 1-5) x 1(T 8 Q. 

This possibility has been explored to some extent in the literature 0, 0, E| , where it has been 
shown that kaon decays place severe constraints on the flavor-changing two-quark couplings of 
the hypothetical new particle, X. It has also been claimed that a light sgoldstino is a viable 
candidate It is well known in the case of light Higgs production in kaon decay that, in addition 
to the two-quark flavor- changing couplings, there are comparable four-quark contributions 0]. They 
arise from the combined effects of the usual SM four-quark \AS\ = 1 operators and the flavor- 
conserving couplings of X. We have recently computed the analogous four-quark contributions to 
light Higgs production in hyperon decay [8j and found that they can also be comparable to the 
two-quark contributions previously discussed in the literature. 

The interplay between the two- and four-quark contributions makes it possible to find models 
with a light Higgs boson responsible for the HyperCP events that has not been observed in kaon 
or 5-meson decay. However, it is not easy to devise such models respecting all the experimental 
constraints. In most models that can generate dsX couplings, the two-quark operators have the 
structure d{l ± 7 5 )sX. Since the part without 7 5 contributes significantly to K — > T\\i r \i~ , their 
data imply that these couplings are too small to account for the HyperCP events 0, In some 
models, there may be parameter space where the four-quark contributions mentioned above and 
the two-quark ones are comparable and cancel sufficiently to lead to suppressed K — > rates 
while yielding S + — > pfi + fi~ rates within the required bounds. However, since in many models the 
flavor- changing two-quark couplings qq'X are related for different (q, q') sets, experimental data on 
5-meson decays, in particular B — > X s /i + /i~, also provide stringent constraints. For these reasons, 
the light (pseudo)scalars in many well-known models, such as the SM and the two-Higgs-doublet 
model, are ruled out as candidates to explain the HyperCP events jsj]. 

In this paper we show that a light pseudoscalar Higgs boson, A\, in the next-to- minimal su- 
persymmetric standard model (NMSSM) |J can be identified with the X particle satisfying the 
constraints 

B{K± -> 7T ± A° l ) < 8.7 x 1(T 9 , B{K S -> vr°A°) < 1.8 x 1(T 9 , 

B(B -> X S A?) < 8.0 x 1(T 7 , (1) 



obtained in R ef. E l from the measurements of K ± 7r ± /x + /i jl^, U|, K$ — > 7i°fi + fi ^|, and 



B — > [13j,|l4j|. We include both two- and four-quark contributions to the K and S + decays, 

neglecting the latter contributions to the B decay. 

In the NMSSM there is a gauge-singlet Higgs field iV in addition to the two Higgs fields H u and 
Hd responsible for the up- and down-type quark masses in the MSSM. We follow the specific model 
described in Ref. with suitable modifications. The superpotential of the model is given by 

W = QY U H U U + QY d H d D + LY e H d E + XH d H u N - ±kN 3 , (2) 
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where Q, U, D, L, and E represent the usual quark and lepton fields, Y u ^ e are the Yukawa couplings, 
and A and k are dimensionless parameters. The soft-supersymmetry-breaking term in the Higgs 
potential is 

K oft = m 2 H jH u \ 2 + m 2 Hd \Hf + m 2 N \N\ 2 - (\A x H d H u N + ±kA k N 3 + K.c.) , (3) 
and the resulting Higgs potential has a global U(1)r symmetry in the limit that the parameters 

There are two physical pseudoscalar bosons in the model which are linear combinations of the 
pseudoscalar components in H u ^ and N. The lighter one, A®, has a mass given by 

m 2 A = 3kxA k + 0(1/ tan/3) (4) 

in the large-tan j3 limit, where x = (N) is the vacuum expectation value of N and tan /3 is the ratio 
of vacuum expectation values of the two Higgs doublets. If the U(1).r symmetry is broken slightly, 
the mass of A® becomes naturally small, with values as low as ~100 MeV phenomenologically 
allowed 0,0]. We will now show that this CP-odd scalar can play the role of the X particle. 

At tree level, the A\ couplings to up- and down-type quarks and to leptons in this model are 
diagonal and described by 

t-Mi = -( l u^uUl 5 u + l d m d d l5 d)^ , C M = l l^I lh lAl (5) 



in the general notation of Ref. [8j used to compute the effect of four-quark operators, where 

/ vS_ l d 

with v = 246 GeV being the electroweak scale and 5_ = (A\ — 2kx)/(A x + kx). Thus, l u can be 
neglected in the large-tan f3 limit. The interactions in C^ qq combined with the operators due to W 
exchange between quarks induce the four-quark contributions discussed earlier. 

For an A\ of mass 214.3 MeV as we propose, the decay into a muon-antimuon pair completely 
dominates over the other kinematically allowed modes: A\ —>■ e + e~, uu, 77. Accordingly, we 
assume B(X —> ~ 1. In addition, the muon anomalous magnetic moment imposes the 

constraint 0] 

\9t\ £ 1-2 • (7) 

This bound results in an A\ width T A < 3.7 x 10~ 7 MeV, which is consistent with the estimate of 
10~ 7 MeV for the width of the Hyper CP p article [f| if \gt\ = \ld\ = 0(1). An Id of order one is also 
consistent with the general estimates [l5| for the size of v5_/ (v2aA- 

At one-loop level, the two-quark flavor- changing Lagrangian for s — > dA\ can be written as 



^Asd = — K d{\ + 75 ) S - m d d(l- l5 )s] A? + H.c. , (8) 



% 9a 

v 

where 

G 



9 A ~ g 7^.2 5^ V qd V qs C Aq 

q=u,c,t 



(9) 
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with Vgr being the elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix and C Aq depending 
on the particles in the loops ^| . The corresponding coupling constant g\ in b — > sA® has a similar 
expression in terms of V* S V b C Aq . 

In the large-tan (3 limit, the pseudoscalar A\ in this model does not couple to up- type quarks and 
Eq. @ is completely determined by parameters in the supersymmetry (SUSY) sector. In particular, 
C Aq can be zero if a super-Glashow-Iliopoulos-Maiani mechanism is operative [151 . In this limit, 
the decays S + — > pA\ and K — > nA® are dominated by four-quark contributions |8| , and the mode 
b — > sA^ is negligible. However, the results of Ref. imply that a vanishing due to C Aq — > 
would not satisfy all the existing constraints. 

This leads us to some details of the model. Assuming that there is only stop mixing, and that 
sup and scharm are degenerate, we have for s — > dA\ or b — > 0], 



C 



5„t tan /3 



-Z d Xj + A m t t> sin 9^ cos 6> f Y u , 



=1,2 



u,u>=l,2 



(10) 



where and are loop-integral functions of parameters in the chargino and squark mixing 
matrices given in Ref. and 9i is the stop mixing angle. Consequently, in this scenario the CKM 
factors involving the top-quark are the relevant ones. It then follows from Eq. that C Aq of the 
right magnitude for £>(£ + — > pA\) ~ 3 x 1CT 8 will yield B(B — ► X S A°) that violates its limit in 
Eq. (JTJ) by about three orders of magnitude. 

Now, the penguin amplitude in this model is dominated by squark and chargino intermediate 
states, and the particular scaling with the top-quark CKM factor in Eq. (jlO[) follows from the 
assumption of stop mixing only. If we assume instead that there is only sup mixing, C A will be 
proportional to 5 qu . In this case, it is possible for A\ to be the Hyper CP particle without violating 
the bounds in Eq. ((TJ). To be specific, assuming that sup mixing is dominant, and that stop and 
scharm are degenerate, we obtain 



C 



Aq 



S qu tan (3 



i=l,2 



«,ui=l,2 



(11) 



In addition, Eq. (jll[) has a vanishing CP-violating phase if all the SUSY parameters are real. This is 
needed if the contributions from two- and four-quark operators are to partially cancel each other 0], 
as needed to satisfy the kaon bounds in Eq. (JTJ). To simplify the discussion, we now consider the 
case where 9 U = 0. In this case 

9a = ^J^i KdVus V2% tan/3j> Xi C/, 2 V a [D 3 (yJ - D 3 (y u>i )] l d , (12) 



where Uij and V%j are the unitary matrices diagonalizing the chargino mass matrix, Dr(v) = 
y \ny/(l — y) is a loop-integral function, and y ri = m~/m?, with m Xi being a chargino mass [15]. 
Therefore, depending on whether the sup mass is larger or smaller than the scharm mass, g A can 
take either sign. 

The expressions for the four-quark contributions to the kaon and hyperon decays are lengthy 
and in general depend on the effective coupling of A\ to a gluon pair 0. For the simple case above, 
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the contributions to A®gg are from quarks in the loop in the large-tan {3 limit, and the squarks 
do not contribute since there is no left- and right-squark mixing. Since L — ► 0, the four-quark 
contributions can be written purely in terms of Id- Using the results in Ref. [8J] (with D — F = 0.25 ), 
the full amplitudes are 



iM{K ± — > 7r ± yl ( j ) ) = g A 



rn 



K 



mz 



(1.08 x 10~ 7 ) l t 



m 



K 



(13) 



PAD 



9 a- 



v -- (6.96 x 10~ 7 ) l d — 



(0.25)^^—^ 



+ (5.72 x 10 6 ) I 1 



v 



ip^ + . (14) 



The corresponding amplitude for K$ 
Ref. ' 



7r A(, which we have not displayed, can be derived from 



Choosing, for example, Id = 0.35, we show the resulting branching ratios in Fig. ^ Even though 
it is possible for A\ to explain the HyperCP observation for a wide range of parameters, as the 
figure indicates, we have found that there are only rather narrow ranges of g A and Id for which the 
kaon bounds are also satisfied, namely 1.0 x 10~ 7 ^ g A /ld ~ 1-3 x 10~ 7 . The 5-meson bound is 
automatically satisfied with our assumptions above. 

We remark that g A in Eq. (fT2J) is proportional to GpV* d V us tan (3m w m Xi , which is of OiA) 
for tan/3 = O(30) and m Xi = C»(100GeV). Thus, in order to have g A = b(10~ 7 ) with these 
parameter choices, one needs f/^Va [D 3 [y ci ) — D 3 (y ui )~\ to be of order 10~ 7 as well. Since as 
mentioned above there is a super-GIM mechanism that can be at work, as long as m„ and mg are 
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FIG. 1: Branching ratios of S + — ► pA\ (solid curve), K + — > -k + A® (dashed curve), and K$ — > 7r A5 
(dotted curve) as functions of g A for Id = 0.35. The horizontal lines indicate the bounds from the HyperCP 
result and the kaon bounds from Eq. (^Q). The allowed range in this case is 0.38 ^ 10 g A 0.43. 
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degenerate enough, one can get an arbitrarily small g A to satisfy the required value. It turns out 
that there is also SUSY parameter space where the desired value of g A can be produced without 
very fine-tuned degeneracy. 

We have carried out a detailed exploration of the parameter space. In general, if the chargino 
masses are larger than the squark masses, fine-tuned squark masses, with splitting of less than one 
GeV between sup and scharm masses, are required. However, if the chargino masses [of (9(100 GeV)] 
are smaller than the squark masses [of (9(1 TeV)], there is some parameter space satisfying the 
constraints with squark-mass splitting of a few GeV to tens of GeV. Taking —Ax = 150 GeV, 
rric = 2.5 TeV, and tan/3 = 30, we display in Fig. EJthe allowed ranges of m s — m s and m 2 /(—Xx) 
corresponding to the required values 1.0 x 10~ 7 g A /ld ~ 1-3 x 10~ 7 , where m 2 is the first element 
of the chargino undiagonalized mass-matrix [15]. In this case, the two chargino eigenmasses are, 
respectively, 100 GeV ^$ m Xl < 150 GeV and 200 GeV ^$ m X2 < 1500 GeV. These chargino masses, 
especially the lighter one, are in the interesting range which may be probed at the LHC and ILC. 
In this figure, we also display the much larger (gray) regions where A\ can explain the HyperCP 
observation. 

In conclusion, we have shown that a light pseudoscalar Higgs in the NMSSM can be identified as 
the possible new particle responsible for the HyperCP events while satisfying all constraints from 
kaon and 5-meson decays. We have found that there are large regions in the parameter space 
where the particle can explain the HyperCP observation. This by itself is not trivial, considering 
that the SM, for example, yields rates that are too large by more than a factor of ten |8j. Only 
after the kaon bounds are imposed does the allowed parameter space become much more restricted. 
We believe that this is sufficiently intriguing to warrant a revisiting of these kaon bounds, perhaps 
by the NA48 experiment. 
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FIG. 2: Parameter space for m^— mg and m 2 /(—Xx) where A\ can explain the HyperCP events (gray 
regions) and simultaneously satisfy the kaon bounds (black regions), with the input described in the text. 
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